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Using the concept of diffused vorticity and the formalism of rotational hydrodynamics we calculate 
the eigenmodes of a harmonically trapped Bose-Einstein condensate containing an array of quantized 
vortices. We predict the occurrence of a new branch of anomalous excitations, analogous to the 
Kelvin modes of the single vortex dynamics. Special attention is devoted to the excitation of the 
anomalous scissors mode. 



PACS numbers: 03.75.Lm 



Introduction 



The existence of the macroscopic wave function de- 
scribing a quantum fluid imposes a velocity flow curl free 
everywhere except on singularity lines known as vortices. 
Around these lines, the circulation of the velocity is non 
zero and is quantized in unit of h/m where m is the 
mass of the particles of the fluid. Recent experiments 
have demonstrated the nucleation of such quantized vor- 
tices in stirred gaseous Bose-Einstein condensates (BEC) 

Multiple quantized vortices are energetically unstable 
in harmonic traps so that for large rotation frequencies, 
Bose-Einstein condensates nucleate several singly quan- 
tized vortices that were observed to form regular triangu- 
lar lattices known as Abrikosov lattices 0,1a, E| . In this 
configuration, the circulation of the velocity field over 
a circle orthogonal to the rotation axis and of radius 
R much larger than the vortex interspacing, is simply 
irR 2 n v h/m, where n v is the density of the vortex lines. 
This is the same formula as that of the velocity field of a 
rigid body rotating at the angular velocity 



^0 = n v — . 
2m 



(1) 



The dynamical properties of a single vortex line were 
first studied by Lord Kelvin H and his results were trans- 
posed to quantum fluids For an excitation of 
wave vector k propagating along the rotation axis, the 
dispersion relation ojk of these modes (Kelvin modes or 
kelvons) is: 



Hk 2 

WK „ i n (i/fc£). 
2m 



(2) 



where £ = (87rpa) -1 / 2 is the healing length giving the 
vortex core diameter, a the scattering length character- 
izing atom binary interactions and p is the density of the 
gas. These modes present the very peculiar feature that 
they only exist with a single helicity as recently proved 



in the experiments of 0| ■ Indeed, the Kelvin-Helmholtz 
theorem that constrains the vorticity to move along with 
the fluid imposes an angular momentum equal to —h to 
the Kelvin modes. 

Although the problem of studying the dynamics of a 
vortex array seems more involved at first sight, it is con- 
siderably simplified for long wavelength perturbations. 
Indeed, in this case a coarse grain averaging method per- 
mits to smooth the discrete nature of vortices. In the 
case of a homogeneous condensate, it was shown that ex- 
citations of wavevector k propagating transversally to the 
rotation axis satisfy the Tkatchenko dispersion relation 
[Hill: 



MIq , 2 



Am 



where Qq is the effective angular frequency of the conden- 
sate defined by The Tkatchenko modes are elastic 
excitations of the lattice and have been recently investi- 
gated theoretically also in the presence of harmonic traps 
[l4|. First experimental evidences for such modes has 
been reported in (lif. 

In this paper we use the coarse grain method to study 
the eigenmodes of a rotating Bose-Einstein condensate 
confined by a harmonic trap using a fully hydrodynamic 
approach including vorticity |16| . In addition to the usual 
collective modes exhibited by the condensate in the ab- 
sence of rotation, we identify an additional branch, anal- 
ogous to the Kelvin modes exhibited by a single vor- 
tex line. With respect to the Tkatchenko modes, whose 
frequencies vanish in the Thomas Fermi limit |14j . the 
Kelvin excitations emerging from our hydrodynamic pic- 
ture approach a finite value in the Thomas-Fermi limit. 
In the case of elongated traps their frequencies actually 
scale like oj 2 /oj±, for a fixed value of Qq/uj±. 

The paper is organized as follows: in Sect. B we de- 
velop the formalism of rotational hydrodynamics in the 
presence of harmonic trapping and derive the general 
equations (|10fl for the dispersion law of the linearized 
excitations. In Sect. C we briefely summarize the results 
for the surface excitations which represent a natural gen- 
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eralization of the modes exhibited by non rotating Bose- 
Enstein condensates. Section D is devoted to the scissors 
mode where an anomalous mode of Kelvin nature is pre- 
dicted. The properties of the scissors oscillations for a 
rotating condensate are discussed in detail using the for- 
malism of linear response theory. Finally in Sect. E we 
obtain a general dispersion relation for the Kelvin modes 
in elongated traps. 



namical equations, including rotational terms, namely: 



d t p = -V(/Ov) 
raftv = -V (ll t + gp + ) — m(V X v) x v, 



A stationary solutions of equations (jSJ) is given by: 



(6) 



B. Rotational hydrodynamics and elementary 
excitations 

It is well known that in the so-called Thomas-Fermi 
regime, where the mean field interaction dominates over 
the quantum pressure, the dynamics of a non rotating 
condensate can be described by the classical equations of 
hydrodynamics: 



dtp = -V(pv) (3) 
m<9 t v = -V (U t + gp + mv 2 /2) , (4) 

where p is the local particle density, g = Aith 2 a/m is the 
coupling constant characterizing the interatomic force 
and v the velocity field satisfying the irrotational con- 
dition V x v = 0. In the case of cylindrical harmonic 
trapping, 



c/ t ^K(.* 2 + y 2 ) + c 2 V) , 

the hydrodynamic equations © and 10} admit a class of 
analytic solutions |17| whose frequencies have been con- 
firmed experimentally with high accuracy. 

In the presence of vortex lines the hydrodynamical for- 
malism must be modified. Here we will employ a simpli- 
fied procedure, by assuming that the characteristic wave- 
length of excitations is large enough so that we can av- 
erage all the physical quantities over domains containing 
several vortices. In this case the average lab frame veloc- 
ity field v is no longer curl-free. Since each vortex carries 
a flux h/m, the average vorticity of the flow is: 



- V x v h 
n = = n v — u, 
2 2m 



(5) 



where n v is the vortex surface density and the unit vec- 
tor u is the local direction of the vortex lines. According 
to equation © , the average vorticity f2 characterizes the 
local vortex distribution: its direction indicates the ori- 
entation of the vortex lines while its modulus is propor- 
tional to the vortex density. 

Like in 0, 0] , we shall simply assume that the av- 
erage velocity field and density satisfy classical hydrody- 



v 

gpo/m 



fl x r 

p - (w 2 



fig) (a*+y*)/2-u>*z*/2, 



where p is the chemical potential and f2o = ^o u z 
characterizes the vortex density and direction accord- 
ing to equation (jSJ). As expected, the stationary ve- 
locity field is equivalent to that of a rigid body ro- 
tating at angular velocity S~2o, while the density is de- 
scribed by the usual Thomas- Fermi profile with the trap- 
ping potential Ut corrected by the centrifugal potential 
U c = —mQ^x 2 + y 2 )/2. The static behavior described 
above has been verified experimentally. In particular, 
the modification of the transverse trapping is used to 
measure experimentally the effective angular velocity f2 
0, @ . More interesting effects concern the study of the 
dynamics of the condensate. In [18j the hydrodynamic 
equations have been used to study the time evolution 
of a condensate containing a vortex array, following the 
sudden switch on of a static deformation in the plane of 
rotation. This produces peculiar non linear effects that 
have been experimentally observed in [y. In the following 
we will focus on the behaviour of the linearized solutions 
which can be derived by looking for small perturbations 
of the the density and velocity field: 



P = Po + Sp 

v = flo x r + <5v, 



with respect to the equilibrium values (in this case, <5v 
can be interpreted as the velocity in the rotating frame). 
If the characteristic wavelength of the perturbation is 
much larger than the vortex spacing, the "averaged" hy- 
drodynamical equations (|5J| can still be applied and, in 
the frame rotating at the angular velocity Jlo , they read, 
in linear approximation, 



d t Sp = V'(p 5v) 

, ( g5p 



d t 5v 



-V 



\ m 



2tt x 5v, 



(8) 
(9) 



where V' denotes the derivation with respect to the co- 
ordinates in the rotating frame. The system of equations 
(|8I9|) constitutes the starting point of our analysis. We 
will study its eigenmodes and show that the spectrum, 
in addition to the usual "phonon" modes displays a new 
class of solutions carrying negative angular momentum 
and that can be physically regarded as Kelvin excita- 
tions. 
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Let us set dt = —ioj' in equations I|8I9[I . Using equation 
(|9*|). we can then express §v as a function of dp as 



<5v = 



1 



1 



/ 2 - 4f2; 



[iw , (f)_ L + 2n xf]- — (f),,, 



where f = —X7'(g5p) jra and (f)|| and (f) ± are, respec- 
tively, the projections of f on the directions parallel and 
orthogonal to z. 

Inserting this expression for the velocity field in the lin- 
earized mass conservation equation © , we get a closed 
equation for the density fluctuations that can be written 



iJ (lo' 2 - 4nl) Sp = A- dp, 
where A denotes the linear operator 



(10) 



A- 8p 



-V, 



— (ioo'VUp + 2Q x V' ± Sp) 
m 



gpo 

m 



,'2 



An 



and po is the equilibrium density profile 0. Equations 
(|1U[) and are the main result of this paper. Simple 
solutions can be found in the form of polynomials. 

Let us stress again that the calculation of lo' is per- 
formed in the rotating frame. The corresponding fre- 
quency lo in the lab frame is obtained through the simple 
relation lo' = lo — m 2 51o, where hm z is the angular mo- 
mentum of the excitation along the rotation axis. Since 
both Sp and its complex conjugate are solutions of ltTT|) 
we find that for any solution with angular momentum m z 
and eigenfrequency lo' , there should be an other solution 



with angular momentum 



and frequency 



Since 



the unperturbed state (condensate rotating at angular 
velocity f2o) is the ground state of the system in the ro- 
tating frame, the physical solutions carrying energy Hlo' 
should correspond to positive frequencies lo' . 

Note also that the equation (|10|l is valid for lo' ^ 
and lo' 2 ^ 4S7q and that these two special cases must 
be treated separately. The case lo' 2 — 417g = does not 
lead to any new mode. The case lo' = is instead more 
interesting since it is related to the Tkatchenko's modes, 
as stressed in the introduction. Starting the analysis back 
from equations © and 10, we can show that the density 
perturbation dp of a zero energy modes depends only 
on the radial coordinate r± — \J x 2 + y 2 and that the 
velocity field is given by: 



2m,n 2 l 



In the case where bp is parabolic, this zero energy mode 
can be simply interpreted as a modification of Qq, i.e. a 
change of the vortex density. 



The rest of the paper is devoted to the solution of 
equation for special cases of physical interest. We 
will first discuss the surface oscillations then the scis- 
sors modes, that will be shown to exhibit a kelvon-like 
branch. Finally we will derive the general dispersion law 
for the Kelvon modes working in the geometry of elon- 
gated traps. 



C. Surface modes 

Surface modes are characterized by the form Sp = 
(x ± iy) 1 and carry angular momentum m z = ±1. In- 
serting this Ansatz in the differential equation l|l(J|) yields 
the following equation for the eigenfrequencies lo' ±1 in the 
rotating frame: 



LO±i ± 2fl Lo' ±l — I (lO a 



o. 



The positive solution of this equation reads: 



J ±l = ^iloI -(i- m 2 t n , 

or, in the laboratory frame, 



ioj\ - (i - i)n 2 ± (i - i)n 



Contrarily to what happens in the case of a non- 
rotating condensate, the two modes are no longer de- 
generate, and the degeneracy lift amounts to: 



ALOl = LO + l - LO-i = 2(1 - 



(11) 



This result is valid for any / [2(| but can be checked for 
some special cases previously reported in the literature. 

In the I = 1 case (dipole motion) one finds Alo = 0. In 
this case the eigenfrequencies are actually unaffected by 
the rotation since lo±i = lo±_ irrespectively of the value 
of fio- This is not surprising since we know that the 
generalized Kohn theorem [lij implies that the center of 
mass motion is determined only by the frequency of the 
harmonic trap. 

The case I — 2 is also well documented, both theo- 
retically and experimentally 0, . Using sum-rule ap- 
proach, it can be shown in particular that the degeneracy 
lift amounts to: 



Alo 



(4 



(12) 



£ z is the angular momentum per particle along z of the 
unperturbed condensate. Since the unperturbed velocity 
field is rigid-like, the angular momentum is given by the 



classical formula (i z 



in (1X21) yields result ITTJ) for I = 2. 



m(r^) fio- Inserting this relation 
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trivial cubic equation: 
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FIG. 



1: Absolute values of the lab frame frequencies of the 
m z — +1 (full line), m z = — 1 (dashed line) and anomalous 
(dash-dotted line) scissors modes for three trap geometries, 
(a): cigar shaped, with u>± = 10cj z . (b): isotropic trap with 
uj z — lo±. In this case the frequency of the anomalous mode 
is zero, (c): pancake trap with uj z = v8u;x. 



D. 



Scissors modes 



The scissors modes are associated with a density fluc- 
tuation 8p — (x ± iy)z carrying angular momentum 
m z = ±1 along the z axis. Using equation one 
finds that the eigenfrequencies uj' must satisfy the non 



± 2fW 2 - w' (uj 2 z +uj 2 ± - Ql) =f 2fi Wz = 0. (13) 



This equation is analogous to that found by the au- 
thors of using the tensor virial method. Two of the 
resulting modes can be identified with the usual m z = ±1 
scissors previously studied in the case of non-rotating 
Bose-Einstein condensates j2^|. The splitting between 
these two modes is easily calculated when fio 0. In 
the laboratory frame one finds Auj = 2VLnU) 2 / {uj 2 + lo\ ) 
in agreement with the sum rule result of [2l| that was 
confirmed experimentally in |25| . This splitting is at the 
origin of the gyroscopic effect investigated in |24j in the 
presence of a single vortex line. The third mode pre- 
dicted by H13|) has instead no analog for firj = 0. A sign 
analysis shows that the frequency uj' a of this anomalous 
mode is positive for m z = — 1 . Just like the kelvons, this 
new scissor mode can only exist with negative helicity. 
The analog of this mode in the case of a single vortex 
line was investigated by |l0j |. The solutions of equation 
(|13|l as a function of f^o are shown in Fig.l for different 
trapping geometries. The interpretation of this mode is 
straightforward for an isotropic trap: in this case, the 
anomalous mode is associated with an overall rotation of 
the condensate and of its lattice. In an isotropic trap this 
rotation indeed costs no energy and the frequency in the 
lab frame is exactly zero. 



For very deformed traps, i.e. for \lo 2 



> ^o, the 



anomalous mode is associated to a scissors mode of the 
vortex lattice while the density profile remains almost at 
rest. 

In the second part of this section we discuss how the 
scissors modes, and in particular the anomalous mode, 
can be excited and observed by suddenly tilting the trap- 
ping potential in the (x, z) plane. To this purpose we 
make use of the formalism of linear response theory. Let 
us apply a perturbing potential of the form 

SU = ez'{x' + iy')e- lJt + e*z'(x' - iy')e iJt 

in the rotating frame. Using the hydrodynamical equa- 
tions, we find that the density and the vorticity will be 
perturbed as 



Sp = e\ u ,z'(x' + iy')e- lul ' t + e*\l,z'{x' -iy')e lul '\ 
5CI = e f i u , 1 e- iu ' t u + + e* l j,Z,e iw ' t u_, 

where u± = Ua; ± iu y . The quantities A and p can be 
calculated using the hydrodynamical equations and one 
finds the result 



/4y 



2nouj 2 z + uj' (uj 2 ± + uj 2 z - nl) 



gP(co>) 



mP{io'Y 



(14) 
(15) 
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where P(w') = uj'^ + 2£Iqu' 2 -uj' (uj 2 z +uj\ -fi§)-2fi w* is 
a polynomial whose roots fix the eigcnfrcqucncies of the 
system (see Ea. (jl3l ). Notice that the change Stto in the 
vorticity, lying in the (x, y) plane, actually corresponds 
to a rotation of the vector lattice. This rotation could 
be imaged experimentally, allowing, together with the 
changes in the shape of the atomic cloud, for a direct 
identification of the various scissors modes. 

On the one hand, the expectation value (z'(x' — iy')) = 
J d rSpz'(x' — iy') can be written as 



{z'(x' - iy')) = ex 
and one obtains the result 



HDg-iw't 



X, 



HD 



— 



(16) 



for the hydro dynamic response function x2 D - 

On the other hand, this expectation value can be calcu- 
lated using the formalism of quantum mechanics. Using 
first order perturbation theory the linear response func- 



tion relative to the operator F = X^kLi z 'k( x 'k ~ *ffe) can 
be in fact written as 



QM 



E 



|(n|Ft| )|2 \(n\F\0)f 



hu' -EL 



hu' + EL 



(17) 



where |0) and |n) are respectively the ground state and 
the excited states of the many body system and E' n are 
the associated eigenenergies in the rotating frame. From 
(|T7|l , we see that the poles of Xu> are the eigenfrequencies 
of the system. Moreover, the sign of the residue at the 
pole gives its helicity: positive residues are associated 
with positive helicity (i.e. with modes excited by z(x + 
iy) ), and negative residues are instead associated with 
negative helicity. 

By identifying and x HD an d by close examination 
of eqs. (|14fl . we can conclude that : 

1. The anomalous mode has positive energy in the 
rotating frame and its angular momentum is m — 
-1. 

2. For lo z < u>±, we have oj a < . For cigar geometries 
the anomalous mode is then associated to a ther- 
modynamical instability of the vortex lattice in the 
laboratory frame [2(j . 

3. On the contrary, for uj z > u>j_, we have cj a > 0. 
For pancake geometries, the vortex lattice is stable 
versus scissor perturbations. 

Let us now consider the special case of the sudden tilt- 
ing at t = of the longitudinal axis of the trap in the lab 
frame. This excitation leads to the following perturbing 
potential: 
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FIG. 2: Density coupling of the scissor m z = +1 (full line), 
m z = —1 (dashed line) and anomalous (dash dotted line) 
transitions for different trapping geometries, (a) Elongated 
trap (u>± = 10cJz) and (b) pancake trap (lo z — \^8ljj_). ) 



5U = axzY (t) = az(x + iy) + az(x — iy)Y(t), 

where Y is the Heavyside step-function equal to for 
t < and 1 otherwise. Expressed in the rotating frame, 
SU reads: 



SU = (az(x + iy)e in ° l + az(x - iy)e~ inot ) Y(t). 
The Fourier transform of this perturbation then yields: 

e u > = a(nS(uj' + fi ) + ij(J + Q )). 

According to equation (|14(l and (|15f) . the induced den- 
sity and vorticity changes are given by: 

Sp = (^j dui'e^X^e-^'*) z(x + iy) + c.c. (18) 
Ml = (^J dLu'e^^e-^'^j u + +c.c. (19) 
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3: Change of vorticity of the scissor m z = +1 (full line), 
m 2 = —1 (dashed line) and anomalous (dash dotted line) 
transitions for different trapping geometries, (a) Elongated 
trap (u>± = WuJz) and (b) pancake trap (u z = v8wj_). For 
these modes the change of vorticity is associated with a change 
of direction of the vortical lattice. ) 



very weak, both in the pancake and cigar traps (except 
for Qq ~ uj±). The behavior of the vorticity is dramati- 
cally different, as observed on Fig. Indeed, while the 
anomalous mode remains weak in the case of an elon- 
gated trap (Fig. 10a)), we see that it dominates the 
dynamics in the pancake geometry. This effect should be 
detectable experimentally by imaging the orientation of 
the vortex lattice. 



E. Kelvin spectrum 

The study of the scissor modes has revealed the exis- 
tence of an anomalous kelvon-like mode with angular mo- 
mentum m z = — 1. For very elongated traps 3> u) z ), 
the corresponding solution is given by: 



2fl 



n 2 



showing that the frequency of the anomalous mode goes 
to zero in the limit of small oj z . In what follows, we 
shall generalize this result by looking for more general 
solutions of (|TT1) satisfying oj' ~ J 1 when ui z — > 0. In 
this approximation, we can restrict equation (flip to non 
vanishing terms in w and ui z , which yields the simplified 
equation 



v ± [2 Po n x v' ± 8 P ] + v: 



* 1 4 3v & 



0. (20) 



Let us now write the density perturbation as the most 
general polynomial of order p + I associated to angular 
momentum m, = —I: 



Sp- 



These integrals can be calculated using standard inte- 
gration techniques in the complex plane. In particular, 
bp and SCI present terms oscillating at the driving fre- 
quency Oo, as well as at the scissor frequencies lo±\ and 
u) a . Let us introduce the Fourier components of these 
modes: 



5p = \Sp e mot + Sp k e- iu '" t ]z(x + iy) + c.c, 

k=±l,a 



(Int 



SSI = <5ft e 



^ Stoke-**'** I lu + cc. 
fc=±l,a 



According to equations l(TS|) and {TSJ, Spk and SClk are 
given, respectively, by the residues of A w < / (u)' + £lo) and 
Plu' /(w' + fio) ■ We have plotted the density response Spk 
on Fig. (|2J|. We see that the anomalous mode weight is 



where the expression is restricted to terms of leading or- 
der in the (rj_, 9, z) cylindrical coordinates and where q m 
is the highest q such that p — 2q > 0. According to equa- 
tion (|2Ufl. the coefficients b q> o must satisfy the recursive 
relation: 



(p-2g + 2)(p+l-2 g )n§ 

9 (-m wV + (p- 2 9)(p- 2 9 + 1 )^) 9_1 ' 

while, for q = 0, 

^-m o u J l+p(p+l)^u J l^bo = 0. (21) 

In order to get non vanishing b q , the coefficient of b 
in equation l|21|l must cancel out. This ensures the quan- 
tization of the eigenfrequencies, since we must have: 



J =p{p+l) 



n§)- 



(22) 
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This dispersion relation represents the generalization 
of the classical Kelvin law to the case of a rotating gas 
confined in a harmonic trap. Like kelvons, these new 
modes have always a negative angular momentum and 
possess a quadratic behavior for large p. If R z and R± 
denote the condensate radii in the longitudinal and trans- 
verse directions respectively, these two quantities are re- 
lated by (see 0) 

u 2 R 2 = (uj 2 ± -n 2 )R 2 ± . 

Moreover, Sp is a polynomial of order p in z. The 
quantity k — p/R z can then be interpreted as the longi- 
tudinal wave vector of the excitation. Using the relation 
= n v h/2m equation (|22|l yields, for large quantum 
numbers p, 

J = N v ^k 2 , 



where N v = n„7ri?5_ is the number of vortices present 
in the condensate. This dispersion law is very similar 
to equation J5J), except for the logarithmic factor that 
vanished through the averaging procedure. 

Note also that equation (1221 is singular for I — 0. This 
is due to the fact that in this case u/ ~ uj z so terms ne- 
glected in (|20|l must be taken into account. In this case, 
there is no decoupling between the phonon and kelvon 
branches, and no simple analytical expression can be ex- 
tracted. 
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